CHIP-FIRING GAMES, POTENTIAL THEORY ON GRAPHS, AND 

SPANNING TREES 

MATTHEW BAKER AND FARBOD SHOKRIEH 

Abstract. Wc study the interplay between chip-firing games and potential theory 
on graphs, characterizing reduced divisors (G-parking functions) on graphs as the 
^D ' solution to an energy (or potential) minimization problem and providing an algorithm 

to efhciently compute reduced divisors. Applications include an "efficient bijective" 
proof of Kirchhoff's matrix-tree theorem and a new algorithm for finding random 
spanning trees. The running time of our algorithms are analyzed using potential 
theory and we show that the bounds thus obtained generalize and improve upon 
several previous results in the literature. We also extend some of these considerations 
to metric graphs. 
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1. Introduction 

Chip-firing games on graphs arise in several different fields of research: in theoretical 
physics they relate to the "abelian sandpile" or "abelian avalanche" models in the 
context of self-organized critical phenomena (jH ETJ EH]); in arithmetic geometry, they 
appear implicitly in the study of component groups of Neron models of Jacobians of 
algebraic curves (|l5l HQl |6]); and in algebraic graph theory they relate to the study of 
flows and cuts in graphs ([3l [121 US])- We recommend the recent survey article [38] for 
a short but more detailed overview of the subject. 

There is a close connection between chip-firing games and potential theory on graphs. 
In this paper, we explore some new aspects of this interplaju. Conceptually, this con- 
nection should not come as a surprise; in both settings the Laplacian operator plays a 
crucial rule. However, in chip-firing games an extra "integrality condition" is imposed; 
in the language of optimization theory, chip-firing games lead to integer programing 
problems whose associated linear programming relaxations can be solved using poten- 
tial theory on graphs. Our potential theory methods allow us to prove some new results 
about chip firing games and to give new proofs and/or generalizations of some known 
results in the subject. We also show that certain "ad-hoc" techniques in the literature 
are naturally explained or unified by our approach. 

Our main potential-theoretic tool is the energy pairing (see §3.31) . which is a canonical 
positive definite bilinear form defined on the set of divisors of degree zerqj on G. This 
pairing can be computed using any generalized inverse of the Laplacian matrix of G. 
The energy pairing can be used to define two functions Sq and bq (for a fixed vertex q) 
on Div(G) which interact in a useful way with chip- firing moves; after firing a set of 
vertices A (not containing the exceptional vertex q) the value of Sq goes down by at 
least the size of the associated cut (see Proposition 14. ip . and the value of bg goes down 
by exactly the size of A (see Proposition 14.41) . 



This paper represents a thorough revision and improvement upon the second author's preprint 
larXiv:0907.4761, which used potential-theoretic methods in only a cursory way. 

A divisor on a finite graph G is an element X^uevCG) '^viv) of the free abelian group Div(G) on the 
set V{G) of vertices of G. The degree of a divisor is the sum of the a„ over all v. 
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Chip-firing moves induce a natural equivalence relation on Div(G) called linear equiv- 
alence of divisors. If we once again fix a vertex q, then a particularly nice set of repre- 
sentatives for linear equivalence classes is given by the q-reduced divisors (see §4.3p . We 
show that the g-reduced divisor equivalent to a given divisor D can be characterized 
as the unique element of \D\q (the set of divisors D' = ^^a^(w) linearly equivalent to 
V for which a^ > whenever v ^ q) minimizing the functional Sg] see Theorem 14. 111 . 
A similar result holds with Sg replaced by bg] see Theorem 14. 13[ Using this result we 
are able to give a new proof of the fact that there is a unique g-reduced divisor in each 
linear equivalence class. 

In order to check whether or not a given divisor is g-reduced, according to the defini- 
tion (given in §4.31) . one needs to check a certain inequality for all subsets of V{G) \ {q}. 
But there is in fact a much more efficient procedure called Dhar's burning algorithm (af- 
ter Dhar [27]); see §5.1[ Using a modification of Dhar's burning algorithm, it is possible 
to obtain an "activity preserving" bijection between g-reduced divisors and spanning 
trees of G; this was originally discovered (using different terminology) by Cori and Le 
Borgne [2S]- In §5.21 we formulate the Cori-Le Borgne algorithm in the language of 
reduced divisors and give a simplified proof of the correctness of the algorithm. 

We then turn to the problem of computing the g-reduced divisor equivalent to a given 
divisor. Dhar's algorithm shows that one can efficiently check whether a given divisor 
is g-reduced; we show in §5.31 (specifically Algorithm Hj) that one can efficiently find 
the unique g-reduced divisor equivalent to a given divisor as well. Algorithm H] can be 
viewed as the search version of Dhar's decision algorithm. The main challenge here is 
the running-time analysis; we use potential theory (specifically, the function bg) to give a 
bound on the running time of Algorithm HI As we have already mentioned, the key point 
is that after firing a set A C V{G) \ {g}, the value of bg goes down by exactly the size 
of A; this makes the function bg a powerful tool for running-time analysis in chip-firing 
processes. The seemingly different techniques of Tardos [IH], Bjorner-Lovasz-Shor [TB] 
or Chung-Ellis [2T], van den Heuvel [21], and Holroyd-Levine-Meszaros-Peres-Propp- 
Wilson [33] all give bounds which are specializations of the running-time bound which 
we derive using bg-, see Remark \5M 

We next turn to applications of Algorithm HI The first application (see §6.11) is an 
"efficient bijective" proof of Kirchhoff's celebrated matrix-tree theorem (stated in a 
more canonical way than usual in Theorem 16. 2p . The efficient bijective matrix-tree 
theorem provides a new approach to the random spanning tree problem, which we state 
in §6.21 This problem has been extensively studied in the literature and there are 
essentially two known types of algorithms for it: determinant based algorithms (e.g. 
[301 [221 [^) and random walk based algorithms (e.g. [HI 111191135] and ^ Chapter 
4]). See Remark 16.31 for possible advantages of our new approach. 

Finally, we study some analagous questions in the context of metric graphs (or "ab- 
stract tropical curves"). While potential theory on finite graphs can be developed purely 
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in the context of linear algebra, potential theory on metric graphs is more conveniently 
formulated in terms of measure theory. We assume in ^ that the reader is familiar with 
some basic facts concerning potential theory on metric graphs as explained in [5l [7]. 
Our main new result is Theorem 17.71 (the metric graph analogue of Theorem I4.1ip : 
using it, we give a new proof of the existence and uniqueness of reduced divisors in 
the metric graph setting. We conclude with a discussion of Dhar's algorithm for metric 
graphs and Luo's search version of this theorem |41j . 



The paper is structured as follows. In ^ we fix our notation. In ^ we recall some 
basic facts from potential theory on graphs and define the energy pairing. The func- 
tionals £q and hq are introduced in §1] and the interplay between chip-firing dynamics 
and potential theory on graphs is studied. The algorithmic applications of this inter- 
play, most notably Algorithm HJ are discussed in §3 Some applications of Algorithm HI 
including an "efficient bijective" proof of Kirchhoff 's matrix-tree theorem and a new 
algorithm for the random spanning tree problem, are discussed in ^ Finally, in ^ we 
extend some our results to metric graphs. 

2. Notation and Terminology 

Throughout this paper, a graph will mean a finite, connected, unweighted multigraph 
with no loop edges. 

The set of vertices of a graph G is denoted by V{G) and the set of edges by E(G). We 
let n = i^V{G) and m = i^E{G). For A C V{G) and v e A, we denote by outdeg^(t;) 
the number of edges of G having v as one endpoint and whose other endpoint lies in 
V{G)\A. 

Let Div(G) be the free abelian group generated by V{G). An element J2vev(G) ^^viv) G 
Div(G) is called a divisor on G. The coefficient ay of (f) in D is denoted by D{v). For 
D e Div(G'), let deg(D) = Y.v&v(g) ^(^) ^nd let Div°(G) be the subgroup of Div(G) 
consisting of divisors of degree zero. We denote by M.{G) = Hom(V^(G), Z) the abelian 
group consisting of all integer- valued functions on the vertices. For A C V{G), xa ^ 
Ai{G) denotes the {0, l}-valued characteristic function of A. Note that {X{v}}v£V{G) 
generates Ai{G) as an abelian group. 

The Laplacian operator A : Ai{G) -^ Div(G) is given by the formula A(/) = 
J2vev{G)^M){v) where 

e={v,w}eE{G) 

This definition naturally extends to all rational- valued functions on vertices. 

Denote by Prin(G') (the group of principal divisors) the image of the Laplacian op- 
erator A : M{G) ^ Div(G). It is easy to see that Prin(G') C Div°(G) and that both 
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Prin(G) and Div°(G') are free abelian groups of rank n — 1. As a consequence, the 
quotient group 

Jac(G') = Div°(G)/Prin(G) 
is finite. Following [3], Ja.c{G) is called the Jacobian of G. 

Let {vi, . . . , v„} be an ordering of V{G). With respect to this ordering, the Laplacian 
matrix Q associated to G is the nxn matrix Q = (qij), where qu is the degree of vertex 
Vi and —Qij {i 7^ j) is the number of edges connecting Vi and Vj. It is well-known (and 
easy to verify) that Q is symmetric, has rank n — 1, and that the kernel of Q is spanned 
by 1, the all-ones vector (see, e.g., [TTl flTj ). 

The ordering {vi, . . . , Vn} of V{G) induces isomorphisms between the abelian groups 
Div(G), A^(G), and the group of n x 1 column vectors with integer coordinates. Under 
these isomorphisms, the Laplacian operator A : M{G) — )■ Div(G) corresponds to the 
matrix Q thought of as a homomorphism Q : Z" — t- Z". More specifically, if [D] 
denotes the column vector corresponding to D G Div(G) and [/] denotes the column 
vector corresponding to / G A4{G), then [A(/)] = Q[f]. 

3. Potential theory 

3.1. Generalized inverses. A matrix has an inverse only if it is square and has full 
rank. But one can define a "partial inverse" for any matrix. 

Definition. Let A be a matrix (not necessarily square). A matrix L satisfying ALA = 
A is called a generalized inverse of A. 

Every matrix A has at least one generalized inverse. In fact more is true: every 
matrix has a unique Moore-Penrose pseudoinversfQ. 

Let Q be the Laplacian matrix of a (connected) graph G. Since Q has rank n — 1, 
it does not have an inverse in the usual sense. But there are several natural ways to 
obtain generalized inverses for Q. Here are some examples. 

Construction 3.1. Fix an integer 1 <i <n and let Qi be the invertible {n — 1) x (n— 1) 
matrix obtained from Q by deleting i^^ row and i^^ column from Q {Qi is sometimes 
called the reduced Laplacian of G with respect to i). Let L(j) be the nxn matrix 
obtained from Q~^ by inserting a row of all zeros after the {i — 1)^* row and inserting 
a column of all zeros after the {i — 1)*^* column. Then L(j) is a generalized inverse of 
Q. Indeed, one checks that QIv(j) = / + -R(j), where I is the nxn identity matrix and 
-R(j) has all —1 entries in the i^^ row and is zero elsewhere; as R{i)Q = 0, we obtain 
QL(i)Q = Q. 



The Moorc-Penrose pseudoinverse of A is a generalized inverse of A having the following additional 
properties: (i) LAL = L and (ii) AL and LA are both symmetric. See [S] for additional details. 
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Construction 3.2. Let J be the n x n all I's matrix. Then Q + Kj \s nonsingular 
and Q^ = (Q + -J)~^ — -J is a generahzed inverse of Q. In fact, Q^ is the unique 
Moore-Penrose pseudoinverse of Q since one easily verifies that QQ^ = Q^Q = I — ^J 
and Q+QQ+ = Q+ . 

One can also use the matrices L(^i^ from Construction 13. ll to obtain other generalized 
inverses for Q: 

Construction 3.3. Let fi = (/ii,/i2,--- , fJ-n)'^ £ K" satisfy X^"=i /^i = 1- Then one 
easily verifies that L^ = Y17=i f^iL{i) is a generalized inverse for Q. The matrix L^ 
has the additional property that L^/i = c^l for some c^ G M; this follows from the 
calculation 



QL^fi = (J + ^ fiiR{i))fJ' = ft- fi 



i=l 



If J is the all-l's matrix as in Construction I3.2[ then G^ = L^ — c^J is also a 
generalized inverse and has the additional property that G^yU = 0. The special case 
where /ij = 1/n for all i gives the Moore-Penrose pseudoinverse Q^ from the previous 
construction. 

As we will see, generalized inverses are intimately related to potential theoretic con- 
siderations on graphs. 

3.2. The j-function. We can think of a graph G as an electrical network in which 
each edge is a resistor having unit resistance. 

Definition. For p,q,v E V{G), let jq{p, v) denote the electric potential at v if one unit 
of current enters a network at p and exits at q, with q grounded (i.e., zero potential). 

From a more mathematical point of view, jq{p,-) is the unique (rational- valued) 
solution to the Laplace equation A/ = (p) — (g) satisfying /(g) = 0; alternatively, one 
can define jq{p, v) to be the {p, f )-entry of the matrix L(q) in Construction 13. II (see, e.g., 
[2U| [S]; note also that (pf ||g) in [12] is the same as our jq{p, v) up to scaling). 

The following properties of the j-function are proved, for example, in [5]: 

• jq{p,q) = 0. 

• 3q{P,v) =Jg{v,p). 

• <Jq{p,v) <jq{p,p). 

• r{p, q) = jq{p,p) = jp{q, q), where r{p, q) denotes the effective resistance between 
p and q. 

3.3. The energy pairing. Let L be any generalized inverse of the Laplacian matrix 
Q. Then the bilinear form (-,•): Div°(G) x Div°(G) -^ Q defined by 

(3.4) {DuD2) = [DifL[D2] 
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is independent of the choice of L. (Indeed, since there are functions /, G Hom(y(G'), Q) 
such that [Di\ = (^[/i] for i = 1, 2, we have (Di, D2) = [fi]'^Q[f2] and it is easy to check 
that the right-hand side does not depend on the choice of /i, /2.) We call the canonical 
bilinear form (■ , ■) the energy pairing on Div°(G). 

Lemma 3.5. The energy pairing is positive definite. 

Proof. Let B be the incidence matrix of the graph. Then Q = BB^, so if [D] = Q[f] 
we have {D,D) = [ffBB^[f] = \\B^[f]\\l D 

Definition. The energy of a divisor D E Div'^(G) is 

£iD) = {D,D) = [DfL[D]. 

Since the energy pairing is positive definite, S{D) > with equality iff Z^ = 0. 

Remark 3.6. The name "energy pairing" comes from the fact that if D G Div'^(G) 
represents an external current in the network, where D{v) units of current enter the 
network at v if D{v) > and —D{v) units of current exit the network at v if D[v) < 0, 
then S{D) is precisely the total energy dissipated in the network. 

We emphasize that the energy pairing is well-defined (i.e., independent of the choice 
of the generalized inverse L) only if the divisors are of degree zero. One can often 
ensure that all divisors under consideration have degree zero by adjusting the value of 
the divisor at some distinguished vertex q. This motivates the definition of the q-energy 
pairing as 

{D,E)^ ={D- deg{D){q),E- deg{E){q)) 

ioi D,Ee'Div{G). 

3.4. The maximum principle. The following result is a graph-theoretic version of 
the classical maximum principle for harmonic functions: 

Lemma 3.7. Let f G A^(G'). Let Amax (resp. A^i^) be the set of vertices where f 
achieve its maximum (resp. minimum) value. Then 

(a) Forv G A^ax, A^(/) > outdeg^^^^ (1;) . 

(b) Forv G Amin, A„(/) < -outdeg^^.^(t'). 

Proof. For part (a) let v G .Amax- For an edge e = vw, if w G Amax then /(t>) = f{w), 
and if w ^ A^ax then f{v) - f{w) > 1. Since A„(/) = T.e=vw&E{G)ifi'^) - /(^))' the 
result follows. Part (b) follows from part (a) by replacing / with — /. 

n 

Corollary 3.8 (Maximum principle). Suppose f G Ai{G) is nonconstant. Then f 
achieves its maximum (resp. minimum) value at a vertex v for which A^(/) > (resp. 
A.(/) < 0;. 

Proof. V can be any boundary point of A^^^ (resp. A^^^^) in Lemma 1X71 D 
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4. Chip-firing dynamics and potential theory 

4.1. Chip- firing dynamics on graphs. Following p], we define an equivalence rela- 
tion ~ (called linear equivalence) on the group Div(G*) as follows: 

Definition. For Di,D2 G Div(G'), Di ~ D2 if and only if Di — D2 is in the image of 
A:M{G) ^Div(G). 

This equivalence relation is closely related to notion of chip-firing games or dollar 
games (see, e.g., [161 Il2l [ISl El [271 H]). Given a divisor D G Div(G), one can view the 
integer D{v) as the number of dollars assigned to the vertex v. 

If D{v) < then v is said to be in debt. A chip-firing move consists of choosing a 
vertex and having it either borrow one dollar from each of its neighbors or give ( "fire" ) 
one dollar to each of its neighbors. For Di,D2 G Div(G), Di ~ D2 if and only if starting 
from the configuration Di one can reach the configuration D2 through a sequence of 
chip-firing moves. 

4.2. Chip-firing moves and the energy pairing. For D G Div(G') we define £q{D) = 
{D, D)q. The following two propositions relate the energy pairing and chip-firing moves, 
and will be used in the next section. 

Proposition 4.1. (a) If E = D + A(/) G Div(G') for some f G M{G), then 

S,{E) = S,{D) + J2 iD + E){v) ■ f{v) -2deg{E) ■ f{q) . 

v€V{G) 

(b) If E = D- A(xa) G Div(G) for some A C V{G)\{q}, then 

S,{E) = S,{D)-J2iD + E){v). 

If, moreover, E is effective on A (i.e. E{v) > for v G A), then 

(4.2) £,{E) < £,{D) - \{A) 

where \{A) denotes the size of the (A, G\A)-cut (i.e., the number of edges having 
one end in A and the other end in V{G)\A). 

Proof, (a) Let deg(i?) = d. Then deg(i3') = d as well. We can find a Q-valued function 
g so that [E + D - 2d{q)] = Q[g]. Then 

£,{E) = {D + A(/) - d{q), D + A(/) - d{q)) 

= S,{D) + {D + E- 2d(g), A(/)) = S,{D) + {Q[9]fLQ[f] 

= £,{D) + [gYQLQ[f] = 8,{D) + [gfQ[f] 

= S,{D) + {Q[g]Y[f] = S,{D) + [D + E- 2d{q)f[f] 

= S,{D) + J2 iD + E){v) ■ f{v) - 2d ■ f{q) . 
vev{G) 
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(b) For the first statement let / = —Xa in part (a). E{v) > for v E A means tliat 
D{v) > A,{xa) = outdeg^(w) for v e A. So E.ga(^ + E){v) > E.eA outdeg^(i;) = 
X{A). D 

Definition. Let 1 denote tlie all-l's divisor. For D G Div(G) and q G V{G), we define 
bg{D) = {l,D)g. If D G Div°(G) we simply write b{D) = {1,D). 

Remark 4.3. bq{D) is tlie "total potential" induced by the external current source given 
by divisor D - deg{D){q) G Div°(G). 

Proposition 4.4. (a) If E = D + A(/) G Div(G') for some f e M{G), then 

(4.5) \{E) = h,{D) + Y, /(^) - (^ - l)/(?) > 

where n = \V{G)\. 
(b) IfE = D- A{xa) e Div(G) for some A C V{G)\{q}, then 

(4.6) b,iE) = \{D) - \A\ , 

where \A\ is the cardinality of the set A. Thus bq{-) is a monovarianl|j. 

Proof, (a) We can find a Q-valued function g so that [1 — n{q)] = Q[g]. Then 
(1,E),= (1,D), + (1,A(/)), 

= (1, D), + iQ[g])^LQ[f] = (1, D), + [gfQLQ[f] 

= {l,D), + [gfQ[f] = {l,D), + {Q[g]nf] 

= {l,D), + [l-n{q)f[f] 

= {l,D), + J2fiv)-{n-l)f{q). 

(b) follows immediately from part (a) for / = —Xa- CD 

The Q-valued function g : V{G) — )■ Q in the proof of Proposition 14.41 (a) can be 
computed explicitly, and this gives a useful formula for bq-. 

Lemma 4.7. Let g^ : V{G) -^ Q be the unique function such that A{gq) = '^^iv)—n{q) 
and gq{q) = 0. Then: 

(a) gq{v) = Y.pev{G)hij>.v) . 

(b) For any divisor D G Div(G), 

(4.8) bq{D) = Y,9MD{v) = ^5^j,(p,^)Z^(^) . 



V p 



In particular, if D{v) > for v ^ q, then bq{D) > 0. 



A quantity which either only goes up or only goes down under some process. 
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Proof. Part (a) is easy and is left as an exercise. Part (b) follows from Construction 13 . ll 
and the definition of the energy pairing. Alternatively, let [1 — n{q)] = Q[gq] and 
deg(L') = d. Then 

{l,D), = {Q[g^]fL[D-diq)] 

= [g,fQL[D - diq)] 

= [g,f[D-d{q)] = J29,iv)D{v) 

V 

V p 

The second statement follows because jq{p, v) > and jq{p, q) = 0. D 

4.3. Reduced divisors. A nice set of representatives for equivalence classes of divisors 
are given by the "reduced divisors" . 

Definition. Fix a vertex q E V{G). A divisor D E Div(G) is called q-reduceaa 
if it satisfies the following two conditions: 

(i) D{v) > for all v E V{G)\{q}. 

(ii) For every non-empty subset A C V{G)\{q}, there exists a vertex v E A such 
that D{v) < outdeg^(t>). 

In other words, every vertex outside q is nonnegative but simultaneously firing all 
the vertices in any non-empty subset A of V{G) which is disjoint from q will result in 
some vertex becoming negative. 

The significance of reduced divisors is mainly due to the fact that for every D E 
Div(G), there is a unique g-reduced divisor D' such that D' ~ D. This basic fact 
was discovered independently (in different guises) by several different authors (see, e.g., 
[28| [26| HU [6]). We give a new proof of this result in Corollary 14. 121 below. 

We wish to study reduced divisors from a potential theoretic point of view. Fix a 
distinguished vertex q and define 

\D\y = {E E Div(G) \E r^ D, E{v) > for all v ^ q} . 

Lemma 4.9. For every D E Div(G) and any vertex q, the set \D\g is non-empty. 

Proof. We can pick an ordering -< on V{G) with the property that q is the first vertex 
in the ordering, and every v ^ q has a neighbor w with w -< v. Now starting from the 
last vertex in this ordering, we can make vertices effective consecutively; at each step 
a vertex v can become effective by replacing D with D — kA{xw) for some neighbor 
w -< V and for some large enough integer k. D 



Reduced divisors are essentially the same thing as G-parking functions |33] or superstable configu- 
rations 1331. 
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Lemma 4.10 (Principle of least action). Let D be a q-reduced divisor. Assume E r^ D 
and write D = E + A(/). 

(a) IfEe \D\^, then f{v) < f{q) for all v G V{G). 

(b) IfE + Aig) e \D\„ then g{q) - g{v) < f{q) - f{v) for all v G V{G). 

Proof. Part (a) is a consequence of Lemma [3171 If / does not achieve its global maximum 
at q, then Amax ^ ^{G)\{q} and for all v G Amax we have A„(/) > outdeg^^^^^^(f ). Since 
D is g- reduced, there must be a vertex u G Amax such that D{u) < outdeg_4_^^^^(M). But 
then E{u) = D{u) — Au{f) < 0, contradicting the assumption that E G \D\q. 

For (b),letE' = E + A(^). Then D = E' + A{f-g) and part (a) gives f{v) - g{v) < 

f{q)-9{q)- □ 

Theorem 4.11. Fix q G V{G) and let D G Div(G). Then D is q-reduced if and only 
ifDe \D\q and£q{D) < Sg{D') for all D' y^ D in \D\g. 

Proof If D is g-reduced, then D G \D\q. Let E G \D\q. Then £^{0) < £q{E); write 
D = E + A{f) with /(g) = 0. By Lemma|4lO](a) we have f{v) < 0. By PropositionO 
(a) we have 

S,{D) = £,{E) + ^(D + E){v) ■ f{v) < S,{E) . 

Now assume D G \D\q and £q{D) < £q{E) for all E G \D\q but D is not g-reduced. 
Then there exists a non-empty set A C V{G)\{q} such that Di = D — A(xa) € \D\q. 
But then by Proposition 14.11 (b) 

8q{Dl) = Sq{D) - X{A) < Sq{D) - 1 . 

It follows that if Sq{Di) = Sq{D2) < Sq{E) ioi all E e \D\q, then both Di and D2 
are g-reduced. By Lemma [4.101 (a) if D2 = Di+ A(/) with /(g) = 0, then f{v) < for 
all V ^ q. Similarly —f{v) > for all f 7^ g, so / = and Di = D2. □ 

Corollary 4.12. Fix q G V{G) and let D G Div(G'). Then there is a unique q-reduced 
divisor D' G Div(G) linearly equivalent to D. 

Proof. It follows from Lemma [4.91 Proposition 14.11 (b), and Lemma [3.51 that we may 
choose D' G Div(G) such that £q{D') < £q{D") for all D" G \D\q. By Theorem [ITU 
D' is the unique g-reduced divisor linearly equivalent io D. D 

An analogue of Theorem 14.111 holds with £q replaced by hq-. 

Theorem 4.13. Fix q G V{G) and let D G Div(G). Then D is q-reduced if and only 
ifDe \D\q and hq{D) < hq{D') for all D' ^ D m \D\q. 

Proof. Note that by Lemma [4.7! (b) the function hq does have a minimum in \D\q. The 
rest of the proof mirrors the proof of Theorem I4.11[ 
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If D is g-reduced, then D G \D\q by definition. Let E G \D\q be anotlier divisor. 
Write D = E + A(/) with /(g) = 0. By Lemma 1110] (a) we have f{v) < 0. Now, by 
Proposition 14.41 (a), we have 

bqiD) = b,iE) + Y,fiv)<ME). 

vjtq 

Now assume D G \D\q and hq{D) < bq{E) for all E G \D\q but D is not g- reduced. 
Then there exists a non-empty set A C V(G)\{g} such that Di = D — A(xa) G l-D|q. 
But then we get a contradiction using the Proposition 14.41 (b) 

hq{Dl) = bq{D) - \A\ < hq{D) - 1 . 

It follows that if hq{Di) = 6g(A) < bq{E) for all E G \D\q, then both D^ and D2 are 
g-reduced. By Lemma Kl (a) if A = ^1 + A(/) with /(g) = 0, then f{v) < for all 
V ^ q. Similarly —f{v) > for all t; 7^ g, so / = and Di = D2. □ 

Remark 4.14. Theorem 14. 131 remains true if bq{D) = {l,D)q is replaced by {h,D)q for 
any "M-divisor" h (i.e. h{v) G M), provided that h{v) > for w 7^ g. 

5. Algorithmic aspects of reduced divisors 

5.1. Dhar's algorithm. Let D be a divisor on the graph G. In order to check whether 
or not D is g-reduced using the definition, one needs to check for all subsets A C 
V{G)\{q} whether or not there is a vertex v E A such that D{v) < outdeg^(u). But 
there is in fact a much more efficient procedure called Dhar's burning algorithm (after 
Dhar 



The idea behind Dhar's algorithm is as follows. Think of the edges of G as being 
made of a flammable material. A flre starts at vertex g and proceeds along each edge 
adjacent to g. At each vertex v 7^ g, there are D[v) flreflghters, each of whom can 
control flres in a single direction (i.e., edge) leading into v. Whenever there are flres 
approaching v in more than D{v) directions, the flre burns through v and proceeds to 
burn along all the other edges incident to v. The divisor D is g-reduced iff the flre 
eventually burns through every vertex of G. 

More formally, Dhar's algorithm can be stated as follows: 

Proof of correctness. If the algorithm outputs FALSE, then D is not g-reduced since 
D{v) > outdeg^^(v) for all v G Ai. If the algorithm outputs TRUE, we claim that 
D is g-reduced; this requires showing that for any subset A C V{G)\{q}, there exists 
a vertex v E A such that D{v) < outdeg^(t'). Let v^ be the minimal element of A 
with respect to the "burning order" on V{G) induced by the algorithm. Then A ^ Ak 
and therefore outdeg^^(wfc) < outdeg^(ffc). But by assumption D{vk) < ontdeg j^^ivk), 
hence D{yk) < outdeg^(ffc) and Vk is the desired element of A. D 
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Input: A divisor D G Div(G), and a vertex q G V{G). 

Output: TRUE if D is g-reduced, and FALSE if D is not g-reduced. 

if D{v) < for some v G V{G)\{q} then output FALSE and Stop. 
Let Aq = V{G) and vq = q. 
for l<z<ra — Ido 

Let Ai = Ai^i \ {vi^i}. 

if for all V G Ai, D{v) > outdeg^.(f) then output FALSE and Stop. 

else let Vi G Ai be any vertex with D{vi) < outdeg^.(t'j). 
end 
Output TRUE. 



Algorithm 1: Dhar's Burning Algorithm 

The complexity of Dhar's algorithm is 0{n'^): there are at most n iterations, and at 
most n inequalities are tested in each iteration. 

5.2. The Cori-Le Borgne algorithm. Using a modification of Dhar's burning al- 
gorithm, it is possible to obtain an "activity preserving" bijection between g-reduced 
divisors (of a given degree d) on G and spanning trees of G. This is more or less just 
a restatement of the work of Cori and Le Borgne in [22] in the language of reduced di- 
visors; however, by using the Cori-Le Borgne algorithm in conjunction with the results 
of ^ and §11 we are able to obtain new results. 

The idea behind the Cori-Le Borgne algorithm is as follows. In our original formu- 
lation of Dhar's algorithm, we burned through multiple edges at once. We now use an 
ordering of E{G) to break ties and implement a "controlled burn" in which only one 
edge at a time is burnt. 

Thus, fix a total order on E{G) and suppose we are given a g-reduced divisor D. We 
run Dhar's burning algorithm on D, starting with a fire at g. However, any time there 
are multiple unburnt edges which are eligible to burn, we always choose the smallest 
one. Whenever the fire burns through a vertex v, we mark the edge along which the 
fire traveled just before burning through v. Since D is g-reduced, the fire eventually 
burns through every vertex of G. The set of marked edges is connected, has cardinality 
n — 1, and covers all vertices and thus forms a spanning tree T^ of G. We thus obtain 
an association {reduced divisors} -^ {spanning trees} (Algorithmic]). 

The remarkable fact discovered by Cori and Le Borgne is that this association is a 
bijection. The inverse map is also completely explicit and can be described as follows 
(using the same total order on E{G)). Suppose we are given a spanning tree T in G. 
A controlled burn starts at the vertex g and, as before, any time there are multiple 
unburnt edges eligible to burn we choose the smallest one. The difference is that 
now the firefighters at v can control incoming fires in every direction except for those 
corresponding to edges of T. Thus the fire burns through a vertex v ^ q exactly when 
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it travels along an edge e e T from some (burnt) vertex w to v. At the moment when v 
is burned through, we set D{v) equal to #{burnt edges adjacent to v} — 1. Eventually 
the fire will burn through every vertex and a nonnegative integer D[v) will have been 
assigned to each vertex v ^ q. The value of D{q) is determined by requiring that 
deg(-D) = d. It turns out that the resulting divisor D is g-reduced, so we obtain an 
association {spanning trees} -^ {reduced divisors} (Algorithm |3]) which one checks is 
inverse to Algorithm [2l 



Input: 

G = (y, E) is graph with a fixed ordering on E, 

q e V{G), 

D = J2v (^v{v), a g- reduced divisor of degree d. 

Output: 

Td a spanning tree of G. 

Initialization: 

X = {q} ("burnt" vertices), 
i? = ("burnt" edges), 
T = ("marked" edges). 

while X ^ V{G) do 

/ = min{e = {s, t} G E{G) \e ^ R, s e X,t ^ X}, 
let V e V{G)\X be the vertex incident to /, 
if a„ = #{e incident to t" | e G R} then 
X^XUiv}, 

T^ru{/}, 

end 

R^RU{f} 
end 

ExtAct= E\R, 
ExtPass= R\T, 
Output To = T. 



Algorithm 2: Reduced divisor to spanning tree. 

Our proof differs from proofs in Cori and Le Borgne; we avoid using their decreasing 
traversals, and we work directly with reduced divisors instead of critical configurations. 

Theorem 5.1. The association given by Algorithms [H and {^ is a bijection. More 
precisely: 

(i) For any q-reduced divisor D of degree d, Algorithm [^ outputs a spanning tree 

Td ofG. 
(ii) For any spanning tree T, Algorithmic outputs a q-reduced divisor Dt of degree 
d on G. 
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Input: 

G = {V, E) is graph with a fixed ordering on E, 

q e V{G), 

T a spanning tree of G. 

Output: 

Dt = J2v^v{v), a g-reduced divisor of degree d. 

Initialization: 

X = {q} ("burnt" vertices), 
i? = ("burnt" edges). 

while X ^ V{G) do 

/ = min{e = {s, t} E E{G) | e ^ i?, s G X, t ^ X}, 
if / G T then 

let V G V{G)\X be the vertex incident to /, 

Ov := #{e incident to w | e G -R}, 

X ^xyj{v} 

end 

R^RU{f} 
end 

ExtAct= E\R, 
ExtPass= R\T, 
Output Dt = ^^CLv{v). 



Algorithm 3: Spanning tree to reduced divisor 

(iii) AlgorithmslE and{3 are inverse to one another: Td^ = T and Dtj, = D. 
Moreover, under the bijection furnished by AlgorithmslE and{^' 
(iv) The set R is the same at the end of both algorithms. 
(v) The externally active edge^ for T are precisely the elements of ExtAct= E\R, 

and the externally passive edges for T are precisely the elements of ExtPass= 

R\T. 
(vi) The degree ofY^^, a^,{v) is equal to g—ex(T), where g = m—n+1. Equivalently, 

Qg = d — g + ex(T). 

Proof, (i) If Algorithm [2] produces an output, it will be a connected subgraph of G with 
n — 1 edges and n vertices, hence a spanning tree. Therefore it suffices to show that if 
X ^ V{G) then {e = {s,t} G E{G) \e ^ R, s e X,t ^ X} is non-empty. For the sake 
of contradiction, assume that {e = {s, t} G E{G) | s G X, t ^ X} C R. Then: 



An edge e G E\T is called externally active for T if it is the largest element in the unique cycle 
contained in T U {e}, and is called externally passive for T if it is not externally active. The external 
activity of T is the number of externally active edges for T and is denoted by ex{T). 
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• Since G is connected, {e = {s, t} e E{G) | s G X, t ^ X} 7^ 0. 

• Let e = {s, t} G E{G) be any edge with s G X and t ^ X. Since s G X, it follows 
from the description of the algorithm that ^{e incident to s | e G R} < a^. 

• The assumption {e = {s, t} G E{G) | s G X, t ^ X} C R implies that outdegjc(s) < 
#{e incident to s\e E R}. 

We conclude that outdegj(^(s) < a^ for any vertex s G dX {dX denotes the set of 
vertices of X that have a neighbor outside X). This contradicts the assumption that 
D is g-reduced. 

(ii) First we show that Algorithm [3] always produces an output (i.e., that it assigns 
a nonnegative integer a„ to each vertex v). Assume X 7^ V{G). Since T is a spanning 
tree of G, there is an edge of T connecting X to V{G)\X. This edge cannot already 
be in R, since otherwise both ends of the edge would have already been added to X. 
This means that Algorithm [3] will not stop until X = V{G), and it always produces an 
output. 

It remains to show that the output is g-reduced. To see this, one can use Dhar's 
algorithm (Algorithm [T]). It is clear that a^ > for all v G V{G)\{q}. Let q = 
vq, vi,- ■ ■ , Vn-1 G V{G) denote the order in which vertices are added to X. Let Aq = 
V{G) and Ai = ylj_i\{wj_i}. Then for any 1 < z < n — 1, we have a„. < outdegJ^.{vi) . 
This is because at least one edge / incident to Vi and connected to V{G)\Ai is added 
to R after a„^ := #{e incident to t'j | e G R} is computed. 

(iii) It suffices to show that Dtj^ = D. Record the sequence of 3-tuples {(X, R^T)} C 
V X E X E which appear as Algorithm [2] runs on the input D and produces the output 
To- Now when Algorithm [3] runs on input Td, the sequence {{X,R,T)} will precisely 
describe the way that sets X and R grow and the edges of the tree To are visited. It 
follows that the value ay which is assigned to each vertex v y^ q in Algorithm [3] matches 
the original input D{v). It follows that a^ and -D(g) are also equal because the degree 
d is fixed. 

(iv) This follows from the proof of (iii) above. Indeed, not only do the sets R match 
at the end of both algorithms, the order in which edges and vertices are "burned" is 
the same in both algorithms. 

(v) When an element / G R\T is added to R, its endpoint v G V{G)\X will be 
revisited and added to X by an element e E T belonging to the fundamental cycle 
associated to / and T. Since e comes after /, we must have f < e (with respect to the 
ordering on E{G)). Therefore / is externally passive. Since / was an arbitrary element 
of R\T, we conclude R\T is precisely the set of externally passive elements of T. 

(vi) A vertex v ^ q is added to X precisely when ay edges incident to v are already 
in R. Subsequently one more edge (corresponding to the edge of the spanning tree) is 
also added to R. In this way, we get a partition of R where each vertex v y^ q accounts 
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for tty + 1 elements of R. As a result, J2v^q i^^ + 1) = #-^- The desired equality 
'^v-iqO'v = 9 — ex{T) follows easily using (v). D 

Remark 5.2. (1) The complexity of both Algorithms [2] and [3] is the same as that of 
Dhar's algorithm (Algorithm [1]) , namely 0{n'^). 
(2) A natural choice for the fixed degree is d = g, in which case it follows from 
Theorem 15.11 (vi) that a^ = ex{T). 

Remark 5.3. The problem of giving an explicit bijection between reduced divisors (in the 
guise of G-parking functions) and spanning trees has been studied in several previous 
works (see, e.g., [251 CHI [10] )• There are also a number of bijections in the literature 
between q-critical configurations and spanning trees (see [281 CSl [IH ISHl HI])- For a 
fixed vertex q, g-critical configurations provide another set of unique representatives 
for equivalence classes of divisors (see, e.g., [121 [IS])- There is a simple relationship 
between reduced and critical divisors: D is g-reduced if and only if K^ — D is g-critical, 
where K+ = T.vgv{g) (deg(iO - l){v) [6]. 

5.3. Computing the reduced divisor. Recall that computing the g-reduced divisor 
equivalent to some divisor D can be viewed as the solution to a linear (Theorem I4.13p 
or quadratic (Theorem 14. lip integer programming problem. Dhar's algorithm (which 
runs in time 0(?t,^)) shows that one can efficiently check whether a given divisor is the 
solution to the corresponding integer programming problem. Next we show that in fact 
one can find the solution efficiently as well. 

Fix a base vertex q G V{G). Given a divisor D G Div(G), Algorithm H] below effi- 
cientlju computes the unique g-reduced divisor D' r^ D. The idea behind the algorithm 
is as follows. Starting with a divisor D, the first step is to replace D with an equiva- 
lent divisor whose coefficients are "small" . This is accomplished by the simple trick of 
replacing [D] with [D] — (5[-^(g)[-D]J, where Lq is as in Construction 13 . 1 1 and [-J denotes 
the coordinate-wise floor function. The second step is to make the divisor effective out- 
side q. This is done by having negative vertices borrow from their neighbors until there 
is no more negative vertex outside q. The third step is to iterate Dhar's algorithm until 
we reach a g-reduced divisor. More specifically, if D is not yet reduced then by running 
Dhar's algorithm on D we obtain a subset A of V{G) \ {g} such that firing all vertices 
in A once yields a divisor D — ^{xa) which is still effective outside g. Replacing D 
hj D — A(xyi) and iterating this procedure, one obtains (after finitely many iterations) 
a g-reduced divisor. Moreover, the number of iterations can be explicitly bounded in 
terms of the j-function ( §3.2p using formula (14. 8p . 

A formal statement of the resulting algorithm appears below (Algorithm |4]) . 

Correctness of Algorithm [4t 

Efficient" in this context means that the running time will be polynomial in ra and n with only 
log(deg(D))-bit computations involved. 
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Input: 

Q is the Laplacian matrix of the graph G, 
D e Div(G), 
q e V{G). 

Output: 

D' r^ D the unique g- reduced divisor equivalent to D. 

(Step 1) 

Find the generahzed inverse L(g) of Q, as in Construction 13 .1] Compute the 

divisor [D'] = [D]- gL^(,)[D] J. 

(Step 2) 

while there exists v ^^ q with D'{v) < do [D'] ^ [D'] + Q[x{v}\- 

(Step 3) 

Let Aq = V{G), Vq = q, and i = 1. 

while i < n — 1 do 

Let Ai = Ai_i \ {vi^i}. 

if there exists Vi G A such that D'ivi) < outdeg^. (wj) then i 4- i + 1. 

else [D'] ^ [D'] - Q[xa^- Reset i = l. 
end 



Algorithm 4: Finding the Reduced Divisor 



Assume for the moment that the algorithm actually terminates and produces an 
output. It is easy to see that the output would be linearly equivalent to D. Also, the 
output would certainly pass Dhar's algorithm and therefore be g- reduced (in fact, as 
discussed above, one can view Step 3 as an iterated Dhar's algorithm). Therefore, for 
the correctness of the algorithm, we only need to show that it always terminates (which 
follows a posteriori from the efficiency analysis below) . 

Efficiency of Algorithm [4l 

Proposition 5.4. // [D'] = [D] - Q[L(^q)[D]\, then \D'{v)\ < deg(i;) for all v ^ q. 
(Note that D'{q) is determined in terms of D{v) for v ^ q because deg{D') = deg{D).) 

Proof. Recall from Construction 13.11 that QL(^q) = / + -R(g), where I is the identity 
matrix and /2(q) has —1 entries in q^^ row and is zero everywhere else. Therefore 
[D] = QL(^g-)[D] + deg{D) ■ e^, where e^ is the column vector which is zero everywhere 
except at position q, which is L Now 



[D'] = [D]-Q[L^,^[D]\ 

= Q{L^,)[D] - [L^,)[D]\) + deg{D) . 
= Qf + deg{D) ■ eg 



eg 
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where f = L(g)[D] — [L(q)[Z}]J is a vector with entries from the interval [0, 1). It is now 
easy to show that the absolute values of the entries of Qf are bounded by the degree 
of the corresponding vertices. D 

Remark 5.5. Computing the generalized inverse L(g) takes time at most 0{n^), where 
u is the exponent for matrix multiplication (currently u = 2.376 [S])- Notice that 
this computation is done only once. The second computation in Step (1) can be done 
using 0{n^) operations (multiplication and addition). For bit complexity, one can check 
that the denominators appearing in the generalized inverse L(q) are annihilated by the 
exponent of the Jacobian group. The exponent is bounded above by the number of 
spanning trees of the graph. If we allow at most c parallel edges then there are at most 
^n-i . ^n-2 spanning trees. Moreover, one can also show that the absolute value of the 
entries of L(g) are bounded above by -Rmax, the maximum effective resistance between 
any two vertices of the graph. Therefore all integers in the algorithm can be represented 
in 0{n ■ log en) bits. 

Now we will use our potential theoretic techniques to bound the number of chip-firing 
moves in Algorithm HJ As we will see, several different bounds in the literature can be 
obtained as special cases or corollaries of our general potential theory bound. 

Proposition 5.6. (a) Let Di be the output of Step 1 of Algorithm^ Then Step 
2 of Algorithm^ terminates in at most bq{K~^ — Di) borrowing moves, where 

^+ = E.ey(G)(deg(tO-l)(^0- 

(b) Let D2 be the output of Step 2 of Algorithm ^ Then Step 3 of Algorithm |^ 
terminates in at most bq{D2) firing moves. 

(c) Algorithm^ terminates in fewer than 

(5.7) ^^'^]q{p,v)deg{v) 

V p 

chip-firing moves. 

Proof, (a) By Proposition 15. 4^ |Di(t')| < deg(t') for all v ^ q. So for any v ^ q with 
Di{v) < 0, only one borrowing is needed to make the vertex positive. Moreover, the 
resulting positive number will clearly be less than the vertex degree again. This fact, 
together with Proposition l5.4[ guarantees that the output of Step 2 satisfies < -Di(f ) < 
deg(?;) for all v ^ q. The result now follows from Lemma [4.71 and Proposition 14.41 (b); 
the value of bq{-) is at least J2v Tl,pJqiPi v)Di{v) on the input of Step 2, and is at most 
^^, ^ ?q(p, f)(deg(w) — 1). Moreover, with each borrowing bq{-) increases by 1. 

(b) This again follows from Lemma [4.71 and Proposition 14.41 (b). Note that D2{v) < 
deg(u), and that no vertex v ^ q will become negative in Step 3. 

(c) This follows from parts (a) and (b) and the inequalities 

bq{K+ -D,)<2 5^ 5^ J,(p, v) deg(tO , 
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V p 

D 

5.4. Comparison to other techniques in the literature. By basic properties of 
the j-function (cf. Section [372]) . we have 

V p V 

By Proposition 15.61 (specifically (15. 7p ). it follows that Algorithm H] terminates in fewer 
than 

(5.8) 3(n — 1) ^ ^ r{v, q) deg(f ) 

V 

chip-firing moves. 

The bound ( 15. Sp can be computed in matrix multiplication time 0{n^) (currently u = 
2.376) because jq{p,v) is simply the (p, w)-entry of the matrix L(q) in Construction 13. ll 

Remark 5.9. There are several ways to bound the expression in (15. 8 p in terms of more 
familiar invariants of the graph. For example: 

(1) Let -Rmax be the maximum effective resistance between vertices of G and let 
^max be the maximum degree of a vertex in G. Then (15. Sp is bounded above by 

3(n - l)/2niax ^ deg(t;) , 

which is, in turn, bounded above by 

3(n- l?R„,..^Ar 



i-niax^^max 



These estimates gives a factor n improvement over the bound for the running 
time of Algorithm H] which could be derived using the technique in [33] by 
Holroyd, Levine, Meszaros, Peres, Propp, and Wilson. 
(2) Using Foster's theorem, one can show that 

r{v,q)<3 Y, {deg{v) + l)-^ 

veV[G) 

(see, e.g. [231 proof of Theorem 6]). So another upper bound for the running 
time of Algorithm H] is 

9{n-l) Y, (deg(i;) + l)-i5^deg(^;). 

This is a good bound when the graph is close (on average) to being regular. 
If one uses the fact that degree of a vertex is at least the edge- connectivity A 
of the graph, one gets a bound of the form 0{v?m/X) for the running time of 
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Algorithm HI This is (up to constant factors) the bound that one can derive 

from the techniques of van den Heuvel [31] . 

(3) Let Ai (called the algebraic connectivity oi G) be the smallest non-zero eigenvalue 

of Q. Then r(p, q) < j- for every p, g G V{G). Indeed, the proof of the Cauchy- 

Schwarz inequality shows that for any positive semidefinite matrix L with largest 

eigenvalue rj, and for all vectors x and y, \x^ Ly\ < ?7||x||2||y||2; if we apply this 

estimate with x = y = (p) — (q), L = Q^ (cf. Construction I3.2p . then rj = 1/Ai, 

and we obtain 

2 
r{p, q) = {{p) - (g), (p) - (g)) < — . 

Ai 

Therefore an upper bound for the running time of Algorithm H] is 

6(n-l) 



^deg 



^1 ^ 

vjtq 

This is the bound that one can derive from the techniques of Bjorner-Lovasz- 
Shor HH or Chung-Ellis [21]. 
(4) By Rayleigh's monotonicity law, we have -Rmax < d{G), where d{G) denotes the 
diameter of G. Equality holds if and only if G is a path. In fact Rmax is much 
smaller than d{G) in a general graph. Another upper bound for the running 
time of Algorithm H] is 

3in-l)diG)Y,degiv). 

This is the bound that one can derive from the techniques of Tardos |48j . 

Remark 5.10. Items (3) and (4) in the previous remark clarify the relationship between 
the seemingly different approaches of Tardos and Bjorner-Lovasz-Shor. 

6. Some applications of the algorithms 

6.1. Bijective matrix-tree theorem. Kirchhoff's celebrated matrix-tree theorem is 
usually formulated as follows. Let G be a (connected) graph. Following the terminology 
from Construction 13.1] fix an integer 1 < i < n and let Qi be the invertible [n — 1) x 
(n — 1) matrix obtained from the Laplacian matrix Q of G by deleting i^^ row and i^^ 
column from Q. 

Theorem 6.1 (Kirchhoff's matrix-tree theorem [36]). The number of spanning trees in 
G is equal to \ det(Qi)|. 

Our aim in this section is to give an "efficient bijective" proof of Theorem 16.11 In 
order to make sense of this goal, it is useful to reformulate Kirchhoff's theorem in a 
more natural way in terms of the Jacobian groupj Jac(G) = Div°(G)/Prin(G), where 



Although wc do not need this here, it is worth mentioning that the energy pairing descends to a 
non-degenerate Q/Z- valued bilinear form on Jac(G); it is called the "monodromy pairing" in |47j . 
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Div°(G) is the subgroup of Div(G) consisting of divisors of degree zero and Prin(G) (the 
group of principal divisors) is the image of the Laplacian operator A : Ai{G) — > Div(G). 
By elementary group theory (e.g. the theory of the Smith normal form), one sees 
that Jac(G') is the torsion partoof the cokernel oi Q : Z,"" ^ Z", and the order of Jac(G') 
is equal to | det(Qj)|. We may thus reformulate Kirchhoff' s theorem in a more canonical 
way as follows: 

Theorem 6.2 (Kirchhoff 's matrix-tree theorem, canonical formulation). The number of 
spanning trees in G is equal to \ Jac(G)|. Moreover, there exists an efficiently computable 
bisection between elements 0/ Jac(G) and spanning trees of G. 

Note that such a bijection cannot be canonical, as that would imply the existence 
of a distinguished spanning tree in G corresponding to the identity element of Jac(G), 
but it is clear (think of the case where G is an n-cycle) that there is in general no 
distinguished spanning tree. Therefore, one needs to make some choices to write down 
a bijection. 

Proof of Theorem \6.2[ By Corollary 14.121 if we fix a vertex q of G, there is a unique 
g-reduced divisor representing each class in Div°(G). In particular, there is an explicit 
bijection between Div°(G) and the set of g-reduced divisors of degree 0. If in addition 
we choose a total order on E{G), then the algorithms in §5.21 show that there is a 
bijection between g-reduced divisors of degree and spanning trees of G. 

We have shown in §5.31 that the unique g-reduced representative for each class in 
Div'^(G)/ Prin(G) can be computed efficiently. And in §5.2l we proved that the bijection 
between g-reduced divisors of degree and spanning trees of G is also efficient. D 

Kirchhoff's matrix-tree theorem is of course a classical result. The main new con- 
tribution here is to observe that the bijections between reduced divisors and spanning 
trees (as described in §5.21 and Remark |5. 31) . in conjunction with Corollary 14. 121 furnish 
a simple bijective proof of Kirchhoff's theorem, and moreover this bijection is efficiently 
computable. 

6.2. Random spanning trees. The random spanning tree problem has been exten- 
sively studied in the literature and there are two known types of algorithms: deter- 
minant based algorithms (e.g. [301 [221 EZ]) and random walk based algorithms (e.g. 
UHl [D Sa ES] and [121 Chapter 4]). 

Here we give a new deterministic polynomial time algorithm for choosing a random 
spanning tree in a graph G. Although the bound we obtain for the running time of 
our algorithm does not beat the current best known running time 0{n'^) of [22], we 
believe that our algebraic method has some advantages. For example, it is trivial that 
the output of our algorithm is a uniformly random spanning tree, whereas in other 



The full cokernel group gives the full "Picard group" Pic(G') = Div(G)/Prin(G'), which is isomor- 
phic to Z © Jac(G) 
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algorithms (e.g. [22]) this fact is non-trivial and requires proof. See Remark 16.31 for 
another advantage. 

The idea behind our algorithm is very simple. Fix a vertex q G V{G) and a total 
ordering of E{G). The first step in the algorithm is to compute a presentation of Jac(G') 
as a direct sum of cyclic groups; this can be done efficiently by computing the Smith 
normal form for Q. Once Jac(G') is presented in this way, it is clear how to select a 
random element. Having done so, one computes the corresponding g-reduced divisor 
and then the corresponding spanning tree. 

This procedure is formalized in Algorithm |5] below. 



Input: A graph G. 

Output: A uniformly random spanning tree of G. 

(1) Fix a vertex q G V{G) and a total ordering of E{G). 

(2) Compute the Smith normal form of the Laplacian matrix of G to find: 

- invariant factors {rii, ...,77,^}, 

- generators {gi, . . . , g^} for Jac(G) (thought of as elements of Div''(G)). 

(3) Choose a random integer < a^ < n^ — 1 for {1 < i < s). 

(4) Compute the divisor D = ^^^^ Ojgs- 

(5) Use Algorithm m to find the unique g-reduced divisor D' equivalent to D. 

(6) Use Algorithm |2] to find the spanning tree corresponding to D'. 



Algorithm 5: Choosing a uniformly random spanning tree 

To our knowledge, the fastest known Smith normal form algorithm (Step (2)) is given 
in [3l] and has running time (tt,^'^^^^^^ log ||(5||)^+°''^\ where ||Q||, for our application, 
means the maximal degree of a vertex Amax- See also |46] for a fast and practical 
Smith normal form algorithm. For the running time of Step (5) see (15.71) . (15. 8p . and 
Remark 15. 91 Step (6) can be done in O(n^) steps. 

Remark 6.3. Note that for repeated sampling of random spanning trees in G, one has 
to perform steps (1) and (2) of Algorithm [5] only once. Note also that if there are A^ 
spanning trees in G, our algorithm uses only log2 N random bits for generating each 
random spanning tree. Thus our algorithm may have some advantages over existing 
methods for sampling multiple spanning trees. For example, very few random bits are 
required in our algorithm to generate pairwise independent spanning trees; to generate 
k pairwise independent spanning trees, the naive approach would use k ■ log2 A^ random 
bits. But one can use standard methods to pick pairwise independent elements of the 
group using only 0(log2 A^) random bits. Also, it is also possible with our method to 
sample multiple spanning trees according to joint distributions other that the uniform 
distribution. (We thank Richard Lipton for these observations, c.f. [39]). 

6.3. Other applications. We list briefiy some other applications of our algorithms: 
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(1) (The group law attached to chip-firing games) If we fix a vertex q G V{G), then 
Jac(G) induces a group structure on the set of g-reduced divisors (G-parking 
functions) or g-critical divisors of G. The latter is called critical group (or 
sandpile group) of G. Adding two elements in one of these groups requires first 
adding the given divisors as elements of Div(G'), and then finding the unique 
g-reduced or g-critical divisor equivalent to the sum. Our algorithm for finding 
g-reduced divisors can be used to efficiently compute the group law in these 
groups. A different approach for performing the group operation is given in 
[3T] using "oil games". The problem of finding a "purely algebraic" method 
for computing the sum of two elements of the critical group (and analyzing 
the running time of the resulting algorithm) was posed as an open problem by 
Chung and Ellis in 



(2) (Determining whether the dollar game is winnable) In [B], the authors consider 
a dollar game played on the vertices of G. Given a divisor D, thought of as 
a configuration of dollars on G, the goal of the game is to get all the vertices 
out of debt via borrowing and lending moves, i.e., to find an effective divisor 
D' linearly equivalent to D. By the proof of Theorem 3.3 in [6], the game 
is winnable iff the unique g-reduced divisor equivalent to D is effective. As a 
corollary, once we can efficiently compute the g-reduced divisor associated to a 
given configuration, we can efficiently decide whether or not there is a winning 
strategy, and when there is one we can efficiently compute a sequence of winning 
moves. 

These considerations are related to the Riemann-Roch theorem for graphs 
from [B]. To any divisor D G Div(G) one associates an integer r{D) > —1, 
called the rank of D, and r{D) > iff the unique g-reduced divisor equivalent 
to D is effective. Our algorithm for finding g-reduced divisors can therefore 
be used to efficiently check whether or not r{D) > 0. More generally, we can 
efficiently check whether r[D) > c for any fixed constant c. It is an open problem 
to determine whether or not one can compute r{D) itself in polynomial time. 



7. Metric graphs 

Our goal in this section is to extend some of the considerations from §3.3[ as well as 
Theorem 14.111 and Corollary 14.121 to the setting of metric graphs. The metric graph 
analogue of Corollary 14.121 has already been proved in |33] and [32], but since this is 
the main ingredient needed to prove the Riemann-Roch theorem for tropical curves (an 
important result in tropical geometry), it seems worthwhile to present the new proof 
which follows. Our main new result is Theorem 17.71 (the metric graph analogue of 
Theorem 14.111) . 
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7.1. Background on metric graphs and potential theory. We assume that the 
reader is famihar with some basic facts concerning potential theory on metric graphs; 
see for example [SI [7]. We recall here the main facts and terminology which we will use. 
Let r be a metric graph. Let C(r) denote the R-algebra of continuous real-valued 
functions on F, let CPA(r) C C(r) be the vector space consisting of all continuous 
piecewise affine functions on F, and let -R(F) be the subgroup of CPA(F) consisting 
of continuous piecewise affine functions with integer slopes (this can be viewed as the 
space of tropical rational functions on F, c.f. [291 USD- 
Let A be the Laplacian operator on F (see [5], [7]), which takes a certain subspace 
BDV(F) of C(F) into the space of measures of total mass zero on F. (The abbreviation 
BDV stands for "bounded differential variation".) More precisely, let Meas'^(F) be the 
vector space of finite signed Borel measures of total mass zero on F and let M C C(F) 
denote the space of constant functions on F. Then the space BDV(F) is characterized 
by the property that / i— )■ A(/) induces an isomorphism of vector spaces 

BDV(F)/M < — > Meas°(F) . 
For a fixed g G F, an inverse to A on {/ G BDV(F) | /(g) = 0} is given by 

(7.1) u^ j j,{x,y)du{y)eBm{r) 

where jq{x,y) is the fundamental potential kernel on F (defined below, see also [7]). It 
is shown in [7] that CPA(F) C BDV(F) and that A(/) is a discrete measure if and only 
if /gCPA(F). 

For / G CPA(F) we have 

A{f) = J2Mf) 
per 

where — o"p(/) is the sum of the slopes of / in all tangent directions emanating from p. 
(Note that A(/) = — div(/) with the conventions from [291132] .) This formula uniquely 
characterizes A on BDV(F), because CPA(F) is dense in C(F) and for / G CPA(F) and 
g G BDV(F), we have 



(7.2) JJA{g) = J^gA{f). 

(In fact, (US) holds for all /,c/ G BDV(F).) 

For fixed q and y, one can define jq{x, y) as the unique element of CPA(F) such that 
Axjqix,y) = ^y{x) — Sq{x). In terms of electrical network theory, if we think of F as 
an electrical network with resistances given by the edge lengths, then jq{x,y) is the 
potential at x when one unit of current enters the network at y and exits at q, with 
reference voltage zero at q. It is a basic fact, proved in [20] (see [5] for an alternate 
proof), that jq{x,y) is jointly continuous in x, y, and q. 
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There is a canonical bilinear form, called the energy pairing, on Meas°(r). It can be 
defined in many equivalent ways; if u = A(f) and u = A(q) are in Meas^fF), then 



It is proved in [3, Theorem 10.4] that the energy pairing is positive definite. If // G 
Meas°(r) is a discrete measure, the energy pairing (yU,yu) has a nice interpretation in 
terms of electrical networks. Write /i = ^ ai6p- — Y^ bj5p. with a^, hj > 0. If a, units 
of current enter the network at each pi and bj units of current exit the network at each 
Qj, then (/i,/i) is the total energy (or power) dissipated in the network. The different 
formulas for the energy pairing in (17. 3p correspond to the different classical formulas 
for computing the power (as the sum of VI, PR, etc.). 

7.2. Existence and uniqueness of reduced divisors. Let Div(r) (the group of 
divisors on F) be the free abelian group on F. We can identify a divisor -D = XI (^i{Pi) 
on F with the discrete measure ^d '■= X]'^«^Pi- We will frequently identify divisors and 
measures and will often not explicitly differentiate between D and fin. 

Let Div°(F) be the subgroup of divisors of degree zero on F, and let Prin(F) be 
the subgroup {A(/) | / G -R(F)} of Div°(F) consisting of principal divisors. We write 
D r^ D' ii D — D' belongs to Prin(F) and say that D and D' are linearly equivalent. 

Fix g G F. A divisor -D = X] ^pip) is called effective if Op > for all p, and is 
called effective outside g if Op > for all p ^ q. For D G Div(F), we define the 
complete linear system \D\ to be the set of all effective divisors E equivalent to D, i.e., 
\D\ = {E G Div(F) \ E > 0,E r^ D}. Similarly, we define \D\q to be the set of divisors 
equivalent to D which are effective outside q: 

\D\g = {Ee Div(F) I E{p) > 0, Vp ^ g, ^ ~ D}. 

Definition. Fix g G F. A divisor Z^ on F is called q-reduced if it satisfies the following 
two conditions: 

(Rl) D is effective outside g. 

(R2) If / G -R(F) is non-constant and has a global maximum at g, then 

D + Aif)^\D\,. 

Note that this is not the usual definition of reduced divisor on a metric graph, but it 
is easily seen to be equivalent to the definition used in [21 |32l Hi] : 

Lemma 7.4. A divisor D is q-reduced if and only if it satisfies (Rl) and 

fR2'j For every closed connected set X C T not containing q, there exists a point 
p G dX such that D{p) < outdegj^ (p). 

Proof. Suppose D is effective outside g and satisfies (R2). Given a closed connected set 
X not containing g, construct a rational function / G -R(F) which is on X and e outside 
of an e-neighborhood of X for some sufficiently small e > 0, with slope 1 in each outgoing 
direction from X and /(g) = e (so / achieves its maximum at g). By (R2), there is 
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a point p e T\{q} such that [D + A{f)){p) < 0. Since (A(/))(p) = -outdeg^b), it 
follows that p G dX and D{p) < outdegx(p)- 

Conversely, suppose D is effective outside q and satisfies (R2'). Given a non-constant 
function / G -R(r) achieving its maximum value at q, let X be a connected component 
of the set of points where / achieves its minimum. By assumption (R2'), there exists 
p G dX such that D{p) < outdegx(p)- Since (A(/))(p) < — outdegx(p), we have 
{D + A(/))(p) < D{p) - outdeg^(p) < and thus D + A(/) ^ \D\g. D 

From the definition, one sees that if D is a g- reduced divisor and D + A(/) is effec- 
tive outside q for some non-constant function / G -R(r), then / cannot have a global 
maximum at q. It turns out that such an / must in fact have a global minimum at q. 
The following lemma is the analogue of Lemma 14.101 

Lemma 7.5 (Principle of least action). If D is q-reduced and f G -R(r) is a rational 
function with D + A(/) G \D\q, then f has a global minimum at q. 

Proof. Suppose not, and let X be the set of points where / achieves its minimum value. 
Then X is a closed connected set not containing q. By Lemma \7A\ there exists p G dX 
such that D{p) < outdegY(p)- On the other hand, we have A(/)(p) < — outdegjY(p), 
and thus {D + A(/))(p) < 0, contradicting the assumption that D + A(/) G \D\q. D 



The importance of reduced divisors is given by the following theorem |32l HS], which 
is analogous to the corresponding result for finite graphs (Corollary 14. 121 above) proved 
in [6]: 

Theorem 7.6. Fix g G F. There is a unique q-reduced divisor in each linear equivalence 
class of divisors on F. 

We will give a new proof of this theorem based on the following energy minimization 
result. For a divisor D G Div(F) of degree d, define the q-energy of D by 

Sq{D) = {D-d{q),D-d{q)) . 

li D = J2^iiPi) is effective and Pi ^ q for all i, then £q{D) is the total energy 
dissipated in an electrical network in which Oj units of current enter the network at Pi 
and d = ^ai units of current exit at q. 

Theorem 7.7. Let D he a divisor on F which is effective outside q. Then D is q-reduced 
if and only if it has smaller q-energy than every other divisor in \D\q. 

Proof. Suppose first that D has minimal g-energy among all divisors in \D\q. We want 
to show that D is g-reduced. Let / G -R(F) be a non-constant rational function with a 
global maximum at q. Without loss of generality, we may assume that /(g) = 0. We 
claim that D' := D + A(/) does not belong to \D\q . Indeed, if D were in \D\q , then 
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we would have 

£,{D) = £,{D') -2{D'- d{q), A(/)) + S,{A{f)) 

> S,{D') - 2f{D') 

> £,{D') 

since f(D') < by our assumptions on / and D' . This contradicts the minimahty of 
Sq{D), and thus D' ^ \D\q as claimed. 

For the other direction, we need to prove that if -D is g-reduced and D' = D + A{f) G 
\D\q with / G -R(r) non-constant, then Sq{D') > Sq{D). As before we may assume that 
f{q) = 0. It follows from Lemma [7.51 that f{p) > for all p G F. We thus have 

Sq{D') = £q{D) + 2{D' - d{q), A(/)) + Sq{A{f)) 

>Sq{D) + 2f{D) 
> ^,{D) 

as desired. D 

Corollary 7.8. For every D G Div(F), there is at most one q-reduced divisor equivalent 
to D. 

In order to deduce Theorem 17.61 from Theorem 17.71 we need some auxiliary results. 

Proposition 7.9. Every divisor D on T is equivalent to a divisor which is effective 
outside q. 

Proof. This is proved in [32] and [13] . Here is simple a variant of the proof in [32] (along 
the lines of a similar argument in [6]). Let S be the set of branch points of F together 
with q and all points in the support of D. If p G F\{g}, define \s{j>) (the "level of p" 
relative to S) to be the minimal number of elements of iSyig} contained in a path from 
q to p. We define \s{q) to be —1. Let k be the maximum value of Xs^v) fo^' P £ T- 
For 2 = 0, ■ ■ ■ , fc let S'j = {p G F I Xsip) = "^l- Fo^' i = 0,- ■ ■ ,k — 1, define Ai to be the 
closure of {j> G F | Xs{p) < i} and define Bi = {p E T \ Xs{p) > i}- We have dAi = Si 
and dBi = Si+i. For i = 0, ■■■, k — 1, let fi he the unique element of -R(F) which is 
on Ai and 1 on i?j . Define Pi = A(/j). By definition. 



b^fHv) - S" a^!Hv) 






where K = outdegg^(t>) > 1 and al, = outdeg^.(w) > 1. It is straightforward to verify 
that we can choose positive integers Cq, Ci, ■ ■ ■ , c^-i inductively, starting from Cfc_i and 
working backwards, in such a way that P := X]i=o ^«-^« satisfies D + P & \D\q. D 

Let Meas^(F) be the space of nonnegative Borel measures of total mass d on F. By 
Alaoglu's theorem, Measi(F) is compact in the weak* topology. 
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Lemma 7.10. Let T be a metric graph. 

(a) ///„ converges to f pointwise on T with fn, f G BDV(r), then A(/„) converges 
to A(/) in the weak* topology on Meas°(r). 

(b) Conversely, fix q & T and suppose that A(/„) converges to A(/) in the weak* 
topology on Meas°(r), with fn, f G BDV(r) and /n(g) = f{q) = for all n. 
Then fn converges uniformly to f onV . 

Proof. (1) Since CPA(r) is dense in C(r), we just need to prove that f gA{fn) — )■ 
f gA{f) for every g G CPA(r). But this is equivalent to f fnA{g) — )> f fA{g), which 
follows from pointwise convergence since A{g) is a discrete measure. 

(2) This follows from the continuity of jq{x,y) together with the identities fnix) = 
J j,{x,y)iAfn)iy) and /(x) = J j,ix,y){Af){y). D 

Let D be an effective divisor of degree d on F. We can topologize \D\ by thinking of 
it as a subset of Meas_,_(r) and giving it the subspace topology. 

Theorem 7.11. The space \D\, with its weak* topology as subset of Meas'^{T), is 
compact. 

Proof. Since Meas_,_(r) is a compact Hausdorff space by Alaoglu's theorem, it suffices 
to show that \D\ is a closed subspace. Suppose /i„ G \D\ converges to /x G Meas^(r). 
We want to show that /i G \D\. Since each yU„ is a discrete measure of total mass 
d with nonnegative integer masses, it follows easily from Alaoglu's theorem and the 
"Portmanteau theorem" (see [H Proof of Proposition 9.5]) that /.i is also a discrete 
measure. Therefore // = /id + A(/) for some / G CPA(r). Since /Xn = /^d + A(/„) 
with fn G -R(r) and /i„ — )■ /i, it suffices to prove that if /„ G -R(r), / G BDV(r), and 
A(/n) converges to A(/), then / G -R(r). We may assume without loss of generality 
that fn{q) = f{q) = for all n, so by Lemma lY. 10( b) /„ converges uniformly to / on P. 
But then every directional derivative of /„ converges to the corresponding directional 
derivative of /. Since the slopes of /„ are integers, it follows that the slopes of / are 
integers as well. Thus / G -R(r) and fi E \D\ as desired. D 

Alternatively, one can prove this using the results of |29] as follows. Gathmann and 
Kerber show that l^l can be given the structure of a compact polyhedral complex by 
thinking of it as a subset of the compact space P^*^^ = T'^/Sd (here Sd denotes the 
symmetric group on d letters). The natural map from \D\ (with the Gathmann-Kerber 
topology) to Meas_,_(P) is a continuous injection whose image is \D\ with its weak* 
topology. As the continuous image of a compact space, |-D| is therefore compact in the 
weak* topology. (Note that since a continuous bijection between compact Hausdorff 
spaces is automatically a homeomorphism, it follows that the Gathmann-Kerber and 
weak* topologies on \D\ are the same.) D 

Remark 7.12. One can also deduce Theorem 17.111 from the results in |43j . 
We now give the promised proof of Theorem I7.6[ 
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Proof of Theore'm\7.6[ By Proposition 17.91 and the definition of reduced divisors, we can 



assume witliout loss of generality that D is effective. We endow \D\ with the weak* 
topology as above. The function £^g : |D| — )• M is continuous since 



^qi.D')= I jq{x,y)dfin'-d(q)ix)dij,D'-d{q){y) 




and jq{x,y) is jointly continuous in x and y. Since \D\ is compact, Sg attains its 
minimum value at some effective divisor Dg G \D\. By Theorem 17. 7[ Dg is g-reduced. 
This proves the existence of a g-reduced divisor equivalent to D. The uniqueness is 
Corollary EE □ 

Remark 7.13. The continuity oi Sg : \D\ -^ M. played an important role in the Proof 
of Theorem 17.61 More generally, the energy pairing is continuous as a function from 
Meas°(r) X Meas°(r) to M. If we fix a reference point q eT, this follows from the joint 
continuity of jg(x, y) as a function of x and y together with the formula 



(/^>^) = / / Jqix,y)dn{x)du{y) . 




Remark 7.14. The continuity of the energy pairing, together with the principle of least 
action (Lemma 17. 5p . can be used to give an a new proof of Amini's theorem [2] that 
for D effective, the map from F to \D\ sending a point g G F to the unique g- reduced 
divisor Dg equivalent to D is continuous. 

We conclude this section with an analogue of Theorem 14.131 for metric graphs. For a 
divisor D G Div(F) of degree d, define its "total potential" hg{D) by 

(7.15) bg{.D) := j jg{x,y)dfiD-diq){x)dy , 

where dy denotes the measure on F whose restriction to each edge is Lebesgue measure. 

Theorem 7.16. Let D be a divisor on F which is effective outside q. Then D is 
q-reduced if and only if bg{D) < hg{D') for every divisor D' ^ D in \D\g. 

The proof is similar to the proofs of Theorem 14.131 and Theorem 17. 7[ so we omit it. 

7.3. Dhar's algorithm for metric graphs and Luo's theorem. In this section, we 
briefly discuss some results from Ye Luo's paper [?T] and their relation to the material 
in the present paper. 

The idea behind Dhar's algorithm for metric graphs [HI Algorithm 2.5] is the same 
as for non- metric graphs. Let F be a metric graph, flx a point g G F, and let D be a 
divisor on F which is effective outside g. Think of the metric graph F as being made 
of a flammable material. At every point p G F with p ^ q, there are D{p) flreflghters, 
each of whom can control flres in a single direction leading into p. A flre starts at g and 
proceeds along each direction emanating from g; whenever there are flres approaching 
p in more than D[p) directions, the flre burns through p and proceeds to burn in all 
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directions emanating from p. (In particular, if e is a segment in F consisting entirely 
of points p with D{p) = 0, then a fire starting at one endpoint of e will burn through 
e unobstructed.) The divisor D is g-reduced iff the fire eventually burns through all of 

r. 

For a formal statement of Dhar's algorithm for metric graphs and a proof of correct- 
ness, see 511 Algorithm 2.5 and Lemma 2.6]. 

In order to algorithmically compute the unique g-reduced divisor equivalent to a 
given divisor D, one can proceed along the lines of Algorithm H] above. Starting with 
a divisor D, the first step is to replace D with an equivalent divisor which is effective 
outside q. This is accomplished by following the proof of Proposition 17. 9[ The second 
step is to iterate Dhar's algorithm until we reach a g-reduced divisor. More specifically, 
if D is not yet reduced then by running Dhar's algorithm on D we obtain a proper 
connected subset Y ^ T oi T which contains q (the "burned portion" of F). Choose a 
connected component X of the complement of Y and (as in the proof of Lemma 17. 4p 
choose e > maximal with respect to the property that there is a rational function 
/ G -R(F) which is on X and e outside of an e-neighborhood of X, with constant 
slope 1 in between. By construction, the divisor D + A{f) is still effective outside q. 
Replacing D by D + A{f) and iterating this procedure, one obtains (after finitely many 
iterations) the unique g- reduced divisor equivalent to D. 

For a formal statement of this algorithm and a proof that it terminates after a finite 
number of iterations, see [H], Algorithm 2.12 and Theorem 2.14]. 

In (I7.15P we defined a metric graph analogue of the functional bg from §4l As in 
Proposition 14.4( b). this functional is a "monovariant" which decreases by an explicit 
function of e when D is replaced hj D + A(/) (with / as above): 

b,{D + A(/)) - b,{D) = J J j,{x,y)A.M)dy 

f{x)A^{jg{x,y))dy 
f{x){6y{x) -5q{x))dy 

U\y) - f\q))dy 

= -iKX)e + ^e^), 

where 1{X) is the total length of X and A(X) is the size of the {X, F\X)-cut. 

It would be interesting to give explicit lower bounds for the e's which can appear 
and thus obtain a running-time analysis of Luo's Algorithm 2.12 similar to our Propo- 
sition [5?6t^c) ; we leave this as an open problem for future research. 





32 matthew baker and farbod shokrieh 

References 

[1] D. J. Aldous. The random walk construction of uniform spanning trees and uniform labelled trees. 

SIAM J. Discrete Math., 3(4):450-465, 1990. 
[2] O. Amini. Reduced divisors and embeddings of tropical curves. 2010. Preprint available at 

arXiv:1007.5364vl. 
[3] R. Bacher, P. de la Harpe, and T. Nagnibcda. The lattice of integral flows and the lattice of 

integral cuts on a finite graph. Bull. Soc. Math. France, 125(2):167-198, 1997. 
[4] P. Bak, C. Tang, and K. Wiesenfeld. Self-organized criticality. Phys. Rev. A (3), 38(1):364~374, 

1988. 
[5] M. Baker and X. Faber. Metrized graphs, Laplacian operators, and electrical networks. In Quan- 
tum graphs and their applications, volume 415 of Contemp. Math., pages 15-33. Amer. Math. 

Soc, Providence, RI, 2006. 
[6] M. Baker and S. Norine. Riemann-Roch and Abel-Jacobi theory on a finite graph. Adv. Math., 

215(2) :766-788, 2007. 
[7] M. Baker and R. Rumely. Harmonic analysis on metrized graphs. Canad. J. Math., 59(2):225-275, 

2007. 
[8] M. Baker and R. Rumely. Potential theory and dynamics on the Berkovich projective line, volume 

159 of Mathematical Surveys and Monographs. American Mathematical Society, Providence, RI, 

2010. 
[9] A. Ben-Israel and T. N. E. Greville. Generalized inverses. CMS Books in Mathematics/Ouvrages 

de Mathematiques de la SMC, 15. Springer- Verlag, New York, second edition, 2003. Theory and 

applications. 
[10] B. A. Benson, D. Chakrabarty, and P. Tetali. G-parking functions, acyclic orientations and span- 
ning trees. Discrete Mathematics, 310(8):1340-1353, 2010. 
[11] N. Biggs. Algebraic graph theory. Cambridge Mathematical Library. Cambridge University Press, 

Cambridge, second edition, 1993. 
[12] N. Biggs. Algebraic potential theory on graphs. Bull. London Math. Soc, 29(6):641-682, 1997. 
[13] N. Biggs. Chip-firing and the critical group of a graph. J. Algebraic Combin., 9(l):25-45, 1999. 
[14] N. Biggs. The Tutte polynomial as a growth function. J. Algebraic Combin., 10(2):115-133, 1999. 
[15] N. Biggs and P. Winkler. Chip-firing and the chromatic polynomial. Report LSE-CDAM-97-03, 

Centre for Discrete and Applicable Mathematics, London School of Economics, London, U.K., 

February 1997. 
[16] A. Bjorner, L. Lovasz, and P. W. Shor. Chip-firing games on graphs. European J. Combin., 

12(4):283-291, 1991. 
[17] B. BoUobas. Modern graph theory, volume 184 of Graduate Texts in Mathematics. Springer- Verlag, 

New York, 1998. 
[18] A. Broder. Generating random spanning trees. Foundations of Computer Science, 30th Annual 

IEEE Symposium on, pages 442-447, 1989. 
[19] D. Chebikin and P. Pylyavskyy. A family of bijections between G-parking functions and spanning 

trees. J. Combin. Theory Ser. A, 110(1):31-41, 2005. 
[20] T. Chinburg and R. Rumely. The capacity pairing. J. Reine Angew. Math., 434:1-44, 1993. 
[21] F. Chung and R. B. Ellis. A chip-firing game and Dirichlet eigenvalues. Discrete Math., 257(2- 

3):341-355, 2002. Kleitman and combinatorics: a celebration (Cambridge, MA, 1999). 
[22] C. J. Colbourn, W. J. Myrvold, and E. Neufeld. Two algorithms for unranking arborescences. J. 

Algorithms, 20(2):268-281, 1996. 
[23] D. Coppersmith, U. Feige, and J. Shearer. Random walks on regular and irregular graphs. SIAM 

J. Discrete Math., 9(2):301-308, 1996. 



CHIP-FIRING GAMES, POTENTIAL THEORY ON GRAPHS, AND SPANNING TREES 33 

[24] D. Coppersmith and S. Winograd. Matrix multiplication via arithmetic progressions. Journal of 

Symbolic Computation, 9(3):251 - 280, 1990. Computational algebraic complexity editorial. 
[25] R. Cori and Y. Le Borgne. The sand-pile model and Tutte polynomials. Adv. in Appl. Math., 

30(l-2):44-52, 2003. Formal power series and algebraic combinatorics (Scottsdale, AZ, 2001). 
[26] R. Cori, D. Rossin, and B. Salvy. Polynomial ideals for sandpiles and their Grobner bases. Theoret. 

Comput. Sci., 276(1-2):1-15, 2002. 
[27] D. Dhar. Self-organized critical state of sandpile automaton models. Phys. Rev. Lett., 64(14):1613- 

1616, Apr 1990. 
[28] A. Gabrielov. Abehan avalanches and Tutte polynomials. Phys. A, 195(l-2):253~274, 1993. 
[29] A. Gathmann and M. Kerber. A Riemann-Roch theorem in tropical geometry. Math. Z., 

259(l):217-230, 2008. 
[30] A. Gucnochc. Random spanning tree. J. Algorithms, 4(3):214-220, 1983. 
[31] J. van den Heuvel. Algorithmic aspects of a chip-firing game. Combin. Probab. Comput., 10(6):505- 

529, 2001. 
[32] J. Hladky, D. Krai, and S. Norine. Rank of divisors on tropical curves. 2007. Preprint available 

at arXiv: 0709. 4485. 
[33] A. E. Holroyd, L. Levine, K. Meszaros, Y. Peres, J. Propp, and D. B. Wilson. Chip-firing and 

rotor- routing on directed graphs. In In and out of equilibrium. 2, volume 60 of Progr. Probab., 

pages 331-364. Birkhauscr, Basel, 2008. 
[34] E. Kaltofen and G. Villard. On the complexity of computing determinants. Comput. Complexity, 

13(3-4):91-130, 2004. 
[35] J. A. Kelner and A. Madry. Faster generation of random spanning trees. In Foundations of Com- 
puter Science, 2009., 50th Annual Symposium on, pages 13-21, Oct 2009. 
[36] G. Kirchhoff. Uber die Auflosung der Gleichungen, auf welchc man bci dcr Untersuchung der 

lincarcn Vcrtcilung galvanischer Strome gcfiihrt wird. Ann. Phys. Chem., (72):497-508, 1847. 
[37] V. G. Kulkarni. Generating random combinatorial objects. /. Algorithms, 11(2):185~207, 1990. 
[38] L. Levine and J. Propp. What is . . . a sandpile? Notices Amer. Math. Soc, 57(8):976-979, 2010. 
[39] R. Lipton. A New Approach To Random Spanning Trees. Blog post (2009), available at 

http: //rjlipton.wordpress . com/2009/07/15/a-new-approach-to-random-spajining-trees/) 
[40] D. Lorenzini. Arithmetical graphs. Math. Ann., 285(3):481-501, 1989. 
[41] Y. Luo. Rank-determining sets of metric graphs. J. Comb. Theory Ser. A, 118:1775-1793, August 

2011. 
[42] R. Lyons and Y. Peres. Probability on Trees and Networks. Cambridge University Press, 2011. In 

preparation. Current version available at http://mypage.iu. edu/^rdlyons/. 
[43] G. Mikhalkin and I. Zharkov. Tropical curves, their Jacobians and theta functions. In Curves and 

abelian varieties, volume 465 of Contemp. Math., pages 203-230. Amer. Math. Soc, Providence, 

RI, 2008. 
[44] A. Postnikov and B. Shapiro. Trees, parking functions, syzygies, and deformations of monomial 

ideals. Trans. Amer. Math. Soc, 356(8):3109-3142 (electronic), 2004. 
[45] M. Raynaud. Specialisation du foncteur de Picard. Inst. Hautes Etudes Sci. Publ. Math., (38):27- 

76, 1970. 
[46] D. Saunders and Z. Wan. Smith normal form of dense integer matrices fast algorithms into practice. 

In Proceedings of the 2004 international symposium on Symbolic and algebraic computation, ISS AC 

'04, pages 274-281, New York, NY, USA, 2004. ACM. 
[47] F. Shokrieh. The monodromy pairing and discrete logarithm on the Jacobian of finite graphs. J. 

Math. CryptoL, 4(l):43-56, 2010. 
[48] G. Tardos. Polynomial bound for a chip firing game on graphs. SI AM J. Discrete Math., 1(3):397- 

398, 1988. 



34 MATTHEW BAKER AND FARBOD SHOKRIEH 

[49] D. B. Wilson. Generating random spanning trees more quickly than the cover time. In Proceedings 
of the Twenty-eighth Annual ACM Symposium on the Theory of Computing (Philadelphia, PA, 
1996), pages 296-303, New York, 1996. ACM. 
E-mail address: mbakerOmath.gatech.edu 

E-mail address: shokrieh@math.gatech.edu 

Georgia Institute of Technology, Atlanta, Georgia 30332-0160, USA 



